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Weak convergence of the empirical copula process is shown to hold under the assumption that 
the first-order partial derivatives of the copula exist and are continuous on certain subsets of 
the unit hypercube. The assumption is non-restrictive in the sense that it is needed anyway to 
ensure that the candidate limiting process exists and has continuous trajectories. In addition, 
resampling methods based on the multiplier central limit theorem, which require consistent 
estimation of the first-order derivatives, continue to be valid. Under certain growth conditions 
on the second-order partial derivatives that allow for explosive behavior near the boundaries, the 
almost sure rate in Stute's representation of the empirical copula process can be recovered. The 
conditions are verified, for instance, in the case of the Gaussian copula with full-rank correlation 
matrix, many Archimedean copulas, and many extreme-value copulas. 
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1. Introduction 

A flexible and versatile way to model dependence is via copulas. A fundamental tool for 
inference is the empirical copula, which basically is equal to the empirical distribution 
function of the sample of multivariate ranks, rescaled to the unit interval. The asymp- 
totic behavior of the empirical copula process was studied in, amongst others, Stute [29], 
Ganssler and Stute [10], Chapter 5, van der Vaart and Wellner [32], page 389, Tsuka- 
hara [30, 31], Fermanian et al. [9], Ghoudi and Remillard [15], and van der Vaart and 
Wellner [33]. Weak convergence is shown typically for copulas that are continuously dif- 
ferentiable on the closed hypercube, and rates of convergence of certain remainder terms 
have been established for copulas that are twice continuously differentiable on the closed 
hypercube. Unfortunately, for many (even most) popular copula families, even the first- 
order partial derivatives of the copula fail to be continuous at some boundary points of 
the hypercube. 



This is an electronic reprint of the original article published by the ISI/BS in Bernoulli, 
2012, Vol. 18, No. 3, 764-782. This reprint differs from the original in pagination and 
typographic detail. 



1350-7265 © 2012 ISI/BS 



2 



Example 1.1 (Tail dependence). Let C be a bivariate copula with first-order 
partial derivatives Ci and C2 and positive lower tail dependence coefficient A = 
liinuioC'(u,M)/w > 0. On the one hand, C'i(m,0) = for all u E [0,1] by the fact that 
C{u, 0) = for all u e [0, 1]. On the other hand, Ci{0, v) lini„_^o C{u, v)/u>X>Q for 
all V G (0, 1]. It follows that Ci cannot be continuous at the point (0, 0); similarly for C2- 
For copulas with a positive upper tail dependence coefficient, the first-order partial deriva- 
tives cannot be continuous at the point (1,1). 

Likewise, for the Gaussian copula with non-zero correlation parameter p, the first-order 
partial derivatives fail to be continuous at the points (0, 0) and (1, 1) if p > and at the 
points (0, 1) and (1, 0) if p < 0; see also Example 5.1 below. As a consequence, the cited re- 
sults on the empirical copula process do not apply to such copulas. This problem has been 
largely ignored in the literature, and unjustified calls to the above results abound. A no- 
table exception is the paper by Omelka, Gijbcls, and Veravcrbeke [22]. On page 3031 of 
that paper, it is claimed that weak convergence of the empirical copula process still holds 
if the first-order partial derivatives are continuous at [0, 1]^ \ {(0,0), (0, 1), (1,0), (1, 1)}. 

It is the aim of this paper to remedy the situation by showing that the earlier cited 
results on the empirical copula process actually do hold under a much less restrictive 
assumption, including indeed many copula families that were hitherto excluded. The 
assumption is non-restrictive in the sense that it is needed anyway to ensure that the 
candidate limiting process exists and has continuous trajectories. The results are stated 
and proved in general dimensions. When specialized to the bivariate case, the condition 
is substantially weaker still than the above-mentioned condition in Omelka, Gijbels, and 
Veraverbeke [22]. 

Let _F be a d-variate cumulative distribution function (c.d.f.) with continuous margins 
Fi,...,Fd and copula C, that is, F{x) = C{Fi{xi), . . . , Fd{xd)) for x e R'^. Let Xi, . . .,X„ 
be independent random vectors with common distribution F^ where Xi = {Xn, . . . ,Xid). 
The empirical copula was defined in Deheuvels [5] as 

Cn{u)=F^{F-,\u,),...,Fj{ud)), ue[0,l]'', (1.1) 

where Fn and Fnj are the empirical joint and marginal cdfs of the sample and where F~^^ 
is the marginal quantile function of the jth coordinate sample; see Section 2 below for 
details. The empirical copula C„ is invariant under monotone increasing transformations 
on the data, so it depends on the data only through the ranks. Indeed, up to a difference 
of order 1/n, the empirical copula can be seen as the empirical c.d.f. of the sample of 
normalized ranks, as, for instance, in Riischendorf [25]. For convenience, the definition 
in equation (1.1) will be employed throughout the paper. 
The empirical copula process is defined by 

C„ = V^iCn - C), (1.2) 

to be seen as a random function on [0, 1]'^. We are essentially interested in the asymp- 
totic distribution of C„ in the space £°°{[0, 1]'^) of bounded functions from [0, l]'^ into M 
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equipped with the topology of uniform convergence. Weak convergence is to be under- 
stood in the sense used in the monograph by van der Vaart and WcUner [32] , in particular 
their Definition 1.3.3. 

Although the empirical copula is itself a rather crude estimator of C, it plays a cru- 
cial rule in more sophisticated inference procedures on C, much in the same way as the 
empirical c.d.f. Fn is a fundamental object for creating and understanding inference pro- 
cedures on F or parameters thereof. For instance, the empirical copula is a basic building 
block when estimating copula densities (Chen and Huang [3], Omelka, Gijbels and Ve- 
raverbeke [22]) or dependence measures and functions (Schmid et al. [27], Genest and 
Segers [14]), for testing for independence (Genest and Remillard [12], Genest, Quessy and 
Rcmillard [11], Kojadinovic and Holmes [17]), for testing for shape constraints (Denuit 
and Scaillet [6], Scaillet [26], Kojadinovic and Yan [18]), for resampling (Remillard and 
Scaillet [24], Biicher and Dette [2]), and so forth. 

After some preliminaries in Section 2, the principal result of the paper is given in Sec- 
tion 3, stating weak convergence of the empirical copula process under the condition that 
for every j G {1, . . . , d}, the jth first-order partial derivative Cj exists and is continuous 
on the set {u £ [0, 1]'': < Uj <!}■ The condition is non-restrictive in the sense that it 
is necessary for the candidate limiting process to exist and have continuous trajectories. 
Moreover, the resampling method based on the multiplier central limit theorem proposed 
in Rcmillard and Scaillet [24] is shown to be valid under the same condition. Section 4 
provides a refinement of the main result: under certain bounds on the second-order partial 
derivatives that allow for explosive behavior near the boundaries, the almost sure error 
bound on the remainder term in Stutc [29] and Tsukahara [31] can be entirely recovered. 
The result hinges on an exponential inequality for a certain oscillation modulus of the 
multivariate empirical process detailed in the Appendix; the inequality is a generalization 
of a similar inequality in Einmahl [7] and was communicated by Hidcatsu Tsukahara. 
Section 5 concludes the paper with a number of examples of copulas that do or do not 
verify certain sets of conditions. 

2. Preliminaries 

Let Xi = {Xii, . . . ,Xid), i € {1,2,...}, be independent random vectors with common 
c.d.f. F whose margins Fi,...,Fd are continuous and whose copula is denoted by C. 
Define Uij = Fj{Xij) for i e {l,...,n} and j £ {l,...,d}. The random vectors Ui = 
{Uii, . . . , Uid) constitute an i.i.d. sample from C . Consider the following empirical distri- 
bution functions: for cc € M'^ and for u G [0, 1]'', 
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Here, order relations on vectors are to be interpreted componentwise, and 1a{x) is equal 
to 1 or according to whether x is an element of A or not. Let Xi;n.j < ■ ■ ■ < Xn-.nj 
and Ui;nj < ■ ■ ■ < Un:n,j be the vectors of ascending order statistics of the jth coordi- 
nate samples Xij , . . . , Xnj and Uij , . . . , Unj , respectively. The marginal quantilc functions 
associated to Fnj and Gnj are 

F-liu,) = M{x e R: F,,j{x) > u,} 

Xk:n,j, ii [k - l)/n<Uj <k/n, 
— oo, if Uj = 0; 

G-}{u,) = miiu e [0, 1]: G,y («) > u,} 

Uk:n,j, if (fc - l)/n < itj < fc/n, 
0, if Uj = 0. 

Some thought shows that Xij < F~^{uj) if and only if Uij < G^j{uj), for all i S 
{!,..., n}, j £ {!,..., d} and Uj £ [0,1]. It follows that the empirical copula in equa- 
tion (1.1) is given by 

In particular, without loss of generality we can work directly with the sample Ui, . . . ,Un 
from C . 

The empirical processes associated to the empirical distribution functions G„ and Gnj 
are given by 

a„(u) = v^(G,i(m) - C(u)), an j {u j) = \/n{Gn J {iij)-iLj), (2.1) 
for u E [0, 1]'' and Uj e [0, 1]. Note that a„j(0) = a„j(l) = almost surely. We have 

a„ ^ a (n — > oo) 

in £°°{[0, l]"^), the arrow denoting weak convergence as in Definition 1.3.3 in van der 
Vaart and Wellner [32]. The limit process a is a G-Brownian bridge, that is, a tight 
Gaussian process, centered and with covariance function 

cov(a(w), a{v)) = C{u Av) ~ C{u)G{v), 

for u,v & [0, 1]''; here uAv= (min(ui,ui), . . . ,min(urf, Vd))- Tightness of the process a and 
continuity of its mean and covariance functions implies the existence of a version of a 
with continuous trajectories. Without loss of generality, we assume henceforth that a is 
such a version. 

For j e {1, . . . let Cj be the jth coordinate vector in R'*. For u e [0, l]'^ such that 
0<Uj <l, let 



Q(.).lim^^^ + '^^)-^(") 



/i->o h 

be the jth first-order partial derivative of G, provided it exists. 
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Condition 2.1. For each j G {1, . . . ,d}, the jth first- order partial derivative Cj exists 
and is continuous on the set Vd,j {u e [0, 1]'': < < 1}. 

Henceforth, assume Condition 2.1 holds. To facilitate notation, we will extend the 
domain of Cj to the whole of [0, l]'* by setting 

lims„p!l<!i±iSLl, ituelO.l]'', u, =0, 

Imisup— if u e [0, 1]'', Uj- = 1. 

In this way, Cj is defined everywhere on [0,1]'', takes values in [0,1] (because \C(u) — 
< l^j ^ ''-'jD' ^1^*^ continuous on the set Vdj, by virtue of Condition 2.1. 

Also note that Cj{u) = as soon as Ui = for some i^ j- 



3. Weak convergence 



In Proposition 3.1, Condition 2.1 is shown to be sufficient for the weak convergence of 
the empirical copula process C„. In contrast to earlier results, Condition 2.1 does not 
require existence or continuity of the partial derivatives on certain boundaries. Although 
the improvement is seemingly small, it dramatically enlarges the set of copulas to which 
it applies; see Section 5. Similarly, the unconditional multiplier central limit theorem for 
the empirical copula process based on estimated first-order partial derivatives continues 
to hold (Proposition 3.2). This result is useful as a justification of certain resampling 
procedures that serve to compute critical values for test statistics based on the empirical 
copula in case of a composite null hypothesis, for instance, in the context of goodness- 
of-fit testing as in Kojadinovic and Yan [18]. 

Assume first that the first-order partial derivatives Cj exist and are continuous 
throughout the closed hypercube [0, 1]'^. For u e [0, 1]'', define 

d 

C{u)^a{u)-Y,C,{u)a,{u,), (3.1) 

where aj(uj) = a(l, . . . , 1, Uj, 1, . . . , 1), the variable Uj appearing at the jth entry. By 
continuity of Cj throughout [0, 1]'', the trajectories of C are continuous. From Fermanian 
et al. [9] and Tsukahara [31] we learn that C„ C as ti — > oo in the space £°°([0, 1]''). 

The structure of the limit process C in equation (3.1) can be understood as follows. 
The first term, a{u), would be there even if the true margins Fj were used rather 
than their empirical counterparts Fnj. The terms —Cj(u)aj(uj) encode the impact of 
not knowing the true quantiles Fj'^(uj) and having to replace them by the empirical 
quantiles F~^{uj). The minus sign comes from the Bahadur-Kiefer result stating that 
\^{G~j{uj) — Uj) is asymptotically undistinguishable from —y/n{Cnj{uj) — Uj); see, for 
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instance, Shorack and Wellner [28], Chapter 15. The partial derivative Cj{u) quantifies 
the sensitivity of C with respect to small deviations in the jth margin. 

Now consider the same process C as in equation (3.1) but under Condition 2.1 and 
with the domain of the partial derivatives extended to [0,1]'' as in equation (2.2). Since 
the trajectories of a are continuous and since aj{0) = aj{l) = for each j € {1, . . . , d}, 
the trajectories of C are continuous, even though Cj may fail to be continuous at points 
u e [0, 1]'', such that Uj G {0, 1}. The process C is the weak limit in ^°°([0, 1]'^) of the 
sequence of processes 

d 

C„(u)=a„(u)-^Q(u)a„,(u,), uelO^lf. (3.2) 

The reason is that the map from £°°{[0,1]'^) into itself that sends a function / to / — 
C'j7rj(/), where {TTj{f)){u) = /(I, . . . , l,Uj, 1, . . . , 1), is linear and bounded. 

Proposition 3.1. // Condition 2.1 holds, then, with C„ as in equation (3.2), 

sup |C„(u) — C„(u)| A (n^oo). 

tie[o,i]'^ 

As a consequence, in ^°°([0, 1]''), 

C„ C (n — > oo). 
Proof. It suffices to show the first statement of the proposition. For u e [0,1]'', put 

Rn{u)^\Cn{u)~Cn{u)l U £ [0, I]''. 

If Uj = for some j £ {1, . . . ,d}, then obviously C„(u) — C„(m) ~ 0, so Rn{u) = as well. 
The vector of marginal empirical quantiles is denoted by 

Vn{u)^{G-l{u,),...,G-l{ua)), u&[Q,lY. (3.3) 

We have 

<C^{u)=^{Cn{u)-C{u)) 

= V^{G„(«„(ii)) - C{vn{u))} + ^{C{v^{u)) - C{u)} (3.4) 
= a„(w„(w)) + ^/n■{C{vn{u)) - C{u)}. 

Since a„ converges weakly in to a C-Brownian bridge a, whose trajec- 

tories are continuous, the sequence (a„)„ is asymptotically uniformly equicontinu- 
ous; see Theorem 1.5.7 and Addendum 1.5.8 in van dcr Vaart and Wellner [32]. As 
suPmj6[o,i] l^nj ("i) ~ %l ~^ almost surely, it follows that 

sup \an{vn{u)) ~ an{u)\^Q (n^oo). 
«e[o,i]'' 
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Fix u e [0, 1]''. Put w{t) = u + t{vn{u) - u} and f{t) = C{w{t)) for t € [0, 1]. If w G 
(0,1]'', then Vniu) G (0,1)'', and therefore w(t) G (0,1)'' for all te (0,1], as weU. By 
Condition 2.1, the function / is continuous on [0,1] and continuously differentiable on 
(0, 1). By the mean value theorem, there exists t* ~ tn{u) G (0, 1) such that /(I) — /(O) = 
/'(t*), yielding 

d 

Vr^{CMu)) - C{u)} = J2(^,{w{t*))V^^iG-}iuJ) - u,). (3.5) 
j=i 

If one or more of the components of u are zero, then the above display remains true as 
well, no matter how t* G (0, 1) is defined, because both sides of the equation are equal to 
zero. In particular, if = for sonic k G {1, . . . , d}, then the kth term on the right-hand 
side vanishes because G,^^(0) = whereas the terms with index j ^ k vanish because the 
fcth component of the vector w(t*) is zero, and thus the first-order partial derivatives Cj 
vanish at this point. 

It is known since Kiefer [16] that 

sup \\/n{G^^^{uj) - Uj) + anj{uj)\A-0 {n~>oo). 

«j6[0,l] 

Since < Cj < 1, we find 



sup 

iie[o,i]'' 



d 

Vri{C{vn{u)) - C{u)} + '^Cj{u + t*{Vn{u) - u})anj{Uj) 



4o 



as n — > oo . It remains to be shown that 



sup Dnj{u)A^O (n — > oo) 

tiG[0,l]<' 



for all j G { 1 , . . . , d} , where 

D,,,{u) = \C,{u + t*{v,,{u)-u})-C,{u)\\a„j{u,)\. (3.6) 

Fix e > and S G (0,1/2). Split the supremum over u G [0,1]^* according to the cases 
Uj G [5,1 — S] on the one hand and uj G [0, (5) U (1 — S, 1] on the other hand. We have 

Pr( sup D„j{u) > ej <Ft( sup Dnj{u)>e] 

+ Pr ( sup D„j (u) > ej. 



Since sup„g[Q_]^]d \vn{u) — m| — >■ almost surely, since Cj is uniformly continuous on 

,e[o,i] 



{u G [0, 1]'': (5/2 <Uj <1 — (5/2}, and since the sequence sup„ grg |a„j(uj)| is bounded 
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in probability, the first probability on the right-hand side of the previous display con- 
verges to zero. As \x — y\ < 1 whenever x,y & [0,1] and since < Cj{'w) < 1 for all 
w G [0,1]'', the second probability on the right-hand side of the previous display is 
bounded by 

Pr( sup \anj{uj)\> e). 

^tij6[0,<5)U(l-<54] ^ 

By the portmanteau lemma, the lim sup of this probability as n — >■ oo is bounded by 

Pr( sup ^^]- 

The process aj being a standard Brownian bridge, the above probability can be made 
smaller than an arbitrarily chosen r/ > by choosing 6 sufficiently small. We find 

limsupPrl sup Dnj{u) > e) < rj. 

n^oo ^tiG[0,l]'' 

As r] was arbitrary, the claim is proven. □ 



An alternative to the direct proof above is to invoke the functional delta method as 
in Fcrmanian et al. [9] . Required then is a generalization of Lemma 2 in the cited paper 
asserting Hadamard differentiability of a certain functional under Condition 2.1. This 
program is carried out for the bivariate case in Biichcr [1], Lemma 2.6. 

For purposes of hypothesis testing or confidence interval construction, resampling 
procedures are often required; see the references in the introduction. In Fermanian et 
al. [9], a bootstrap procedure for the empirical copula process is proposed, whereas in 
Rcmillard and Scaillct [24], a method based on the multiplier central limit theorem is 
employed. Yet another method is proposed in Biichcr and Dettc [2]. In the latter paper, 
the finite-sample properties of all these methods are compared in a simulation study, 
and the multiplier approach by Remillard and Scaillet [24] is found to be best over- 
all. Although the latter approach requires estimation of the first-order partial deriva- 
tives, it remains valid under Condition 2.1, allowing for discontinuities on the bound- 
aries. 

Let ^1, 1^2, ■ • • be an i.i.d. sequence of random variables, independent of the random vec- 
tors Xi,X2, . . . , and with zero mean, unit variance, and such that ^/Pi'dCi \> x)dx < 
oo. Define 

1 " 

""(^) = ^E^»(l{^»i ^ FniM,- ■■,X,d< Fjiu,)} - C„(w)). (3.7) 
In (€°°([0, 1]''))^, we have by Lemma A.l in Remillard and Scaillet [24], 



(a„,a^) (a,a') (ti^oo). 



(3.8) 
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where a' is an independent copy of a. Further, let Cnj{u) be an estimator of Cj{u); for 
instance, apply finite differencing to the empirical copula at a spacing proportional to 
n~^/^ as in Remillard and Scaillet [24]. Define 

d 

C',(w) =a^,(M) -^C„j(u)a^j(uj), (3.9) 
i=i 

where a'„j{uj) = . . . , l,Uj, 1, . . . , 1), the variable Uj appearing at the jth coordinate. 

Proposition 3.2. Assume Condition 2.1. If there exists a constant K such that 
\Cnj{u)\ < K for all n, j,u, and if 

sup |(7„,(u)-Q(m)| (n— >oo) (3.10) 

tl(E[0,l]'':«3G[(54-(5] 

for allSe (0,1/2) and all j e {l,...,d}, then in (£°°([0, 1]''))^ we have 

(C„,C)-(C,C') (n^oo), 
where C is an independent copy of C 
Proof. Recall the process in equation (3.7), and define 

d 

C'^{u)=a'„{u)-Y,C,{u)a'^^{u,), € [0, if . 

i=i 

The difference with the process in equation (3.9) is that the true partial derivatives 
of C are used rather than the estimated ones. By Proposition 3.1 and equation (3.8), we 
have 

(C„,C;)^(C,C') (n^oo) 

in (£°°([0,1]''))2. Moreover, 

d 

It suffices to show that each of the d terms on the right-hand side converges to in 
probability, uniformly in w g [0, 1]''. The argument is similar to the one at the end of the 
proof of Proposition 3.1. Pick 6 G (0, 1/2), and split the supremum according to the cases 
Uj g[S,1 — S] and Uj € [0,6) U (1 — <5, 1]. For the first case, use equation (3.10) together 
with tightness of a'„j. For the second case, use the assumed uniform boundedness of 
the partial derivative estimators and the fact that the limit process aj is a standard 
Brownian bridge, having continuous trajectories and vanishing at and 1. □ 
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4. Almost sure rate 

Recall the empirical copula process C„ in equation (1.2) together with its approxima- 
tion Cn in equation (3.2). If the second-order partial derivatives of C exist and are 
continuous on [0,1]^^, then the original result by Stute [29], proved in detail in Tsuka- 
hara [30], reinforces the first claim of Proposition 3.1 to 

sup |C„(m) - C„(-u)| 
Me[o.i]'' 

(4.1) 

= 0{n "'^/'^(logn)"'^/^(loglogn)^/^) (n ^ od) almost surely. 

For many copulas, however, the second-order partial derivatives explode near certain 
parts of the boundaries. The question then is how this affects the above rate. Recall 
Vd,j ={ue [0, 1]'*: 0<Uj< 1} for j € {1, . . . 

Condition 4^.1. For every i,j (z {1, . . . ,d} , the second-order partial derivative Cij is 
defined and continuous on the set Vd^i H Vd.j , and there exists a constant K > such that 

I (w) I < K min ( —-^ , -—^ ) , u^Vd,^f^ Vd,j . 

Condition 4.1 holds, for instance, for absolutely continuous bivariate Gaussian copulas 
and for bivariate extreme-value copulas whose Pickands dependence functions are twice 
continuously differentiable and satisfy a certain bound; see Section 5. 

Under Condition 4.1, the rate in equation (4.1) can be entirely recovered. The fol- 
lowing proposition has benefited from a suggestion of John H.J. Einmahl leading to an 
improvement of a result in an earlier version of the paper claiming a slightly slower rate. 
Furthermore, part of the proof is an adaptation due to Hideatsu Tsukahara of the end 
of the proof of Theorem 4.1 in Tsukahara [30], upon which the present result is based. 

Proposition 4.2. // Conditions 2.1 and 4-1 are verified, then equation (4-1) holds. 

Proof. Combining equations (3.4) and (3.5) in the proof of Proposition 3.1 yields 

d 

with a„ the ordinary multivariate empirical process in equation (2.1), Vn(u) the vector 
of marginal empirical quantiles inequation (3.3), and w(t*) = u + 1* {vn{u) — u} a certain 
point on the line segment between u and w„(it) with local coordinate t* = tn{u) € (0, 1). 
In view of the definition of C„(u) in equation (3.2), it follows that 

sup |C„(m)-C„(u)| <I„ + II„ + III„, 

«e[o,i]'' 
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where 

L„ = sup |q;„(w„(-u)) - a„(u)|, 

u£[0,l]^ 
d 

II„ V sup I ^/^(G^/ [uj)- uj ) + anj {uj ) | , 

d 

Illn^V] sup Dnj{u), 

with Dnj{u) as defined in equation (3.6). By Kicfcr [16], the term II„ is 0(n~^/"'(logri)^/^ x 
(log log n)^/"') as n — > oo, almost surely. It sufEees to show that the same almost sure rate 
is valid for I„ and 111^, too. 

The term I„. The argument is adapted from the final part of the proof of Theorem 4.1 
in Tsukahara [30] , and its essence was kindly provided by Hideatsu Tsukahara. We have 

In<M„(A„), max sup |G,7/(mj) - 

je{i,...,d}„^g[o,i] 

and Mn{a) is the oscillation modulus of the multivariate empirical process q„ de- 
fined in equation (A.l). We will employ the exponential inequality for Pr{M„(a) > A} 
stated in Proposition A.l, which generalizes Inequality 3.5 in Einmahl [7]. Set a„ = 
n~^''^(loglogn)^/^. By the Chung-Smirnov law of the iterated logarithm for empirical 
distribution functions (see, e.g., Shorack and Wellner [28], page 504), 



limsup — sup |G„j'^(uj) — Ujl = limsup — sup \vj—Gnj{vj)\ 

n^oo CLn ujg[0,l] n->oo On vj€[0,l] 

(4.2) 

= l/v2 almost surely. 

Choose A„ = 2/1^2 ^^^n~^/'*(logn)^/^(loglogn)^/'* for K2 as in Proposition A.l. Since 
^n/ {n^^'^CLn) — >■ as n — > 00, and since the function if) in equation (A. 2) below is de- 
creasing with V'(O) = 1; it follows that ■i/;(A„/(ri^/^a„)) > 1/2 for sufficiently large n. 
Furthermore, we have 



E- 



n>2 



2a„ / 4:1 n3/2(loglogn)i/2 

n>2 \ ^ ^ / 

By the Borel-Cantelli lemma and Proposition (A.l), as 71 — ?> 00, 

I„ < Mn{An) < Af„(a„) = 0(n^^^^{\ogn)^/^{\og\ogn)^/^) almost surely. 
The term III„ . Let 

5n = n"^/^(logn)(loglogn)"^/^. 
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Fix j G {l,...,d}. We split the supremum of D„j(u) over u € [0, 1]*^ according to the 
cases Uj G [0, (5„) U (1 — (5„, 1] and Uj G [(5„, 1 — 5„] . 

Since <Cj < 1, the supremum over u G [0, 1]'' such that uj G [0,(5„) U (1 — (5„, 1] is 
bounded by 

sup D„j{u) < sup \anj{Uj)\. 

iie[0,l]'':"je[0,<5„)U(l-(5„,l] Uje[0,(5„)u(l-i5„,l] 

By Theorem 2.(iii) in Einmahl and Mason [8] apphed to (d, z^, fc„) = (l,l/2,n5„), the 
previous supremum is of the order 

sup |a„,(«,)| = 0(<5y2(loglogn)i/2) 

UjG[0,5„)U(l-i5„l 

(4.3) 

= 0(7i-i/''(logn)i/2(iQgiog„)i/4^ ^ 

almost surely. 

Next let u G [0, l]'* be such that (5„ < Uj < 1 — By Lemma 4.3 below and by convexity 
of the function (0, 1) 9 s i-> l/{s(l — s)}, 

D„j{u) = \Cj{u + Xn{u){Vn{u) - u}) - C.j{u)\\anjiUj)\ 

\uj{l- Uj ) G„/ (uj ) ( 1 - G„/ (wj ) ) y 

with |la;||i = l^^j l- Let 5„ = (logn)!/^ jog^Qg^. clearly Y.'^=2'^~^K'^ < °°- ^y 

Csaki [4] or Mason [20], 

Pr I sup , |""^ ''!!Lo > infinitely often ) = 0. 
Vo<s<i (s(l - s))^/^ / 

It follows that, with probability one, for all sufficiently large n, 

\an]{Uj)\ < {Uj{l - Uj))^^^bn, Uj G [0, 1]. 

Let / denote the identity function on [0,1], and let |j-||oo denote the supremum norm. 

For Uj G [Sn, 1 - Sn], 



i-g;/K)>(i-.,)(i- "^"^-^/"°° ). 

By the law of the iterated logarithm (see (4.2)) 

j|G^j^ ^ ^11 = o((5„) [n—^oo) almost surely. 
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We find that with probabihty one, for all sufficiently large n and for all u G [0, l]'^ such 
that Uj G [5„, 1 — 5n], 

D„,{u) < 2K{u,{l - u,))-^^^\\vn{u) - u\\ibn. 

We use again the law of the iterated logarithm in (4.2) to bound ||vn(u) — As a con- 
sequence, with probability one, 

sup Dnj{u) = 0{S~^/^\og\ogn)^/^n-^/^b^) 

lie[0,l]<':M,G[i5„,l-i5„] 

(4.4) 

= 0(n ^/"'(loglogn)^/'*) (n oo) almost surely. 

The bound in (4.4) is dominated by the one in (4.3). The latter therefore yields the total 
rate. □ 

Lemma 4.3. If Conditions 2.1 and 4-1 hold, then 

\C,{v) - C,{u)\ < Xmaxf — i / ) \\v - u^, (4.5) 

for every j G {1, . . . , d} and for every u,v £ [0, l]"* such that < Uj < 1 and < Vj < 1; 
here \\x\\i = l^d denotes the Li-norm. 

Proof. Fix j e {l,...,d} and u,v E [0,1]'' such that Uj,Vj e (0,1). Consider the line 
segment w{t) ~ u + t{v — u) for t e [0, 1], connecting w{0) = u with w{l) = v; put Wi{t) = 
Ui + t{vi — Ui) for i £ {1, . . . ,d}. Clearly < Wj{t) < 1 for all t £ [0, 1]. Next, consider 
the function f{t) = Cj{w{t)) for t G [0,1]. The function / is continuous on [0,1] and 
continuously differcntiable on (0, 1). Indeed, if Ui ^ Vi for some i £ {1, . . . , d}, then < 
Wi{t) < 1 for all t £ (0, 1); if Ui ~ Vi, then Wi{t) ^ Ui ~ Vi docs not depend on t at all. 
Hence, the derivative of / in i G (0, 1) is given by 



iei 



where X ^ {i £ {1, . . . ,d}: Ui ^ Vi}. By the mean- value theorem, we obtain that for some 

t*e(0,i), 

C,iv) - C,{u) = /(I) - /(O) = fit*) = Y,{v. - u.da,iw{t*)). 

iei 

As a consequence, 



\Cj{u) — Cj{v)\ <\\v — u\\iiRax sup \Cij{w{t))\. 

o<i<l 
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By Condition 4.1, 

1 



\Cj{u)-Cjiv)\<\\v~u\\iK sup 



o<t<i Wj{t){l - Wj{t)}' 

Finally, since the function s i— )■ l/{s(l — s)} is convex on (0, 1) and since Wj{t) is a convex 
combination of Uj and Vj, the supremum of l/[wj{t){l — Wjlt)}] over i G [0, 1] must be 
attained at one of the endpoints Uj or Vj. Equation (4.5) follows. □ 



5. Examples 

Example 5.1 (Gaussian copula). Let C be the d-variate Gaussian copula with cor- 
relation matrix R G M'*'^'', that is, 

C{u) = Pr l^fl < u,}^ , 7. e [0, 1]^ 

where X = [Xi, . . . ,Xd) follows a d-variate Gaussian distribution with zero means, unit 
variances, and correlation matrix R; here $ is the standard normal c.d.f. It can be checked 
readily that if the correlation matrix R is of full rank, then Condition 2.1 is verified, and 
Propositions 3.1 and 3.2 apply. 

Still, if < pij ~ corj:{Xi,Xj) < 1 for all j G {2,...,d}, then on the one hand we 
have lim„j^4,oC'i(Mi,u_i) = 1 for all u_i G (0, l]''"^, whereas on the other hand we have 
Ci{u) = as soon as Uj = for some j G {2, . . . ,d}. As a consequence, Ci cannot be 
extended continuously to the set {0} x ([0, 1]'^-^ \ (0, 1]''"^). 

In the bivariate case. Condition 4.1 can be verified by direct calculation, provided the 
correlation parameter p satisfies |p| < 1. 

Example 5.2 (Archimedean copulas). Let C be a c?-variate Archimedean copula; 
that is, 

C{u) = 4>-\4>{ui) + ■■■ + 4>{ud)), u G [0, 1]^ 

where the generator 0: [0, 1] — > [0,cxj] is convex, decreasing, finite on (0, 1], and vanishes 
at 1. whereas 0"-'^ : [0, oo) [0, 1] is its generalized inverse, <t)~^{x) = ini{u G [0,1]: (p{u) < 
x}] in fact, if d > 3, more conditions on (j) are required for C to be a copula; see McNeil 
and Neslehova [21]. 

Suppose (p is continuously differentiable on (0,1] and (f>'{0+) = — oo. Then the first- 
order partial derivatives of C are given by 

^'^""^ ^ /{cilr ^e[o,i]^o<^,<i. 

If Ui = for some i 7^ j, then C{u) = and (f)'{C{u)) = —oo, so indeed Cj{u) = 0. We find 
that Condition 2.1 is verified, so Propositions 3.1 and 3.2 apply. 
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In contrast, Cj may easily fail to be continuous at some boundary points. For instance, 
if 0'(1) = 0, then Cj cannot be extended continuously at (1, . . . , 1). Or if (j>~^ is long-tailed, 
that is, if linix^oo (p'^ix + y)/<p~^{x) = 1 for all y G K, then lim^jj^o C(ui, = 1 

for all U-i € (0, l]'^^^, whereas Ci{u) = as soon as Uj = for some j e {2, . . . ,d}; it 
follows that Ci cannot be extended continuously to the set {0} x ([0, l]''"^ \ (0, 1]''"^). 

Example 5.3 (Extreme-value copulas). Let C be a d-variate extreme-value copula; 
that is, 

C{u) = cxp(-£(- logMi, . . . , - logud)), u e (0, 1]'', 
where the tail dependence function ^■. [Q,ooY [0,oo) verifies 

l{x) = I max {wjXj)H{dw), a;G[0,oo)'', 

with H a non-negative Borel measure (called spectral measure) on the unit simplex 
Ad-i = {w G [0, l]*^: wi + ■ ■ ■ + Wd = 1} satisfying the d constraints / WjH{dw) = 1 for 
all J € {1, . . . , d}; sec, for instance, Leadbettcr and Rootzcn [19] or Pickands [23]. It can 
be verified that £ is convex, is homogeneous of order 1, and that max(a;i , . . . , Xd) < i{x) < 

Xi-\ + Xd for aU x € [0, oo)'*. 

Suppose that the following holds: 



For every j G {1, . . . , d}, the first-order partial derivative £j of £ with respect (5-1) 
to Xj exists and is continuous on the set {x G [0, oo)'^: Xj > 0}. 

Then the first-order partial derivative of C in u with respect to Uj exists and is continuous 
on the set {u G [0, 1]'': <Uj < 1}. Indeed, for u G [0, 1]"^ such that <Uj < 1, we have 

n( N / ~^^j(~log'"i'---'~log"rf)' if > for ah i, 

^0, if = for some i ^ j. 

The properties of £ imply that < £j < 1 for all j G {1, . . . , d}. Therefore, if Ui for 
some i j, then Cj{u) — > 0, as required. Hence if (5.1) is verified. Condition 2.1 is verified 
as well and Propositions 3.1 and 3.2 apply. 

Let us consider the bivariate case in somewhat more detail. The function A: [0, 1] — > 
[1/2, 1] : t H> A{t) = £{1 — t, t) is called the Pickands dependence function of C. It is convex 
and satisfies max(t, 1 — t) < A{t) < 1 for all t G [0, 1]. By homogeneity of the function £, we 
have £{x,y) — {x + y)A{^^) for {x,y) G [0,oo)^ \ {(0,0)}. If A is continuously differen- 
tiable on (0, 1), then (5.1) holds, and Condition 2.1 is verified. Nevertheless, if ^(1/2) < 1, 
which is always true except in case of independence (A = 1), the upper tail dependence 
coefficient 2{1 — ^(1/2)} is positive so that the first-order partial derivatives fail to be 
continuous at the point (1,1); see Example 1.1. One can also see that Ci will not admit 
a continuous extension in the neighborhood of the point (0,0) in case ^'(0) = —1. 
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We will now verify Condition 4.1 under the following additional assumption: 



The function A is twice continuously differentiable on (5 2) 

(0, 1) and M = supo<t<i{i(l - t)A"{t)} < oo. 

In combination with Proposition 4.2, this will justify the use of the Stute-Tsukahara 
almost sure rate (4.1) in the proof of Theorem 3.2 in Genest and Segers [13]; in particular, 
see their equation (B.3). Note that the weight function t(l — t) in the supremum in (5.2) 
is not unimportant: for the Gumbel extreme- value copula having dependence function 
A{t) = {t^/'^ + (1 - ty/^Y with parameter 9 e (0, 1], it holds that A"{t) ^- oo as i ^ or 
i — >■ 1 provided 1/2 < 6* < 1, whereas condition (5.2) is verified for all 6 € (0, 1]. 
The copula density at the point {u,v) G (0, 1)^ is given by 

■• , , Ciu,v) f , , , , in-t)A"it) 



uv \ log(uw) 
where 

t = e (0, 1), ^lit) = A{t) - tA'it), v{t) = A{t) + (1 - t)A\t). 

Note that if A"{l/2) > 0, then Ci2{w, w) — > oo as w 1 1- The properties of A imply < 
fj.(t) < 1 and < v{t) < 1. From - log(a;) > l-x, it follows that -l/log(uu) < min{l/(l- 
u), 1/(1 — u)} for {u,v) £ (0,1)^. Since C{u, v) < min(w, v) and since min(a, b) min(c, d) < 
min{(ac), (bd)} for positive numbers a,b,c,d, we find 

^ / N min(M, v) ( , , . / 1 
< Ci2iu,v) < ^^<^ 1 + i\/min' 

uv 



1-li' 1 



< (1 + i\/)min 
Similarly, for {u,v) G (0, 1) x [0, 1], 



1 1 



u(l — u) ' v{l — v) 



C(u,v) f , , iHl-t)A"{t)\ 

^''1 fi{t){l-^i{t))+ \ . , ' I ifO<z;<l, 



Cu{u,v) = <^ V? \ log(u 

, 0, if V e {0, 1}. 

Continuity at the boundary v = follows from the fact that C{u,v) — )■ as f 0; con- 
tinuity at the boundary v — 1 follows from the fact that t and ii{t) — )■ as v — )• 1. 
Since — log(M) < {1 — u)/u, we find, as required, 

0<-Cn{u,v)<^-^-^^^, (u,«)e(0,l)x [0,1]. 
u[l — u) 

Example 5-4 (If everything fails. . . ). Sometimes, even Condition 2.1 does not hold: 
think, for instance, of the Frechet lower and upper bounds, C(m, v) = max(u + v — 1,0) 
and C{u,v) = miii{u^v), and of the checkerboard copula with Lebesgue density c = 
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21[o,i/2]2u[i/2,i]2 • In these cases, the candidate limiting process C has discontinuous tra- 
jectories, and the empirical copula process does not converge weakly in the topology of 
uniform convergence. 

One may then wonder if weak convergence of the empirical copula process still holds in, 
for instance, a Skorohod-type topology on the space of cadlag functions on [0, 1]^. Such 
a result would be useful to derive, for instance, the asymptotic distribution of certain 
functionals of the empirical copula process, for example, suprcma or integrals such as 
appearing in certain test statistics. 



Appendix: Multivariate oscillation modulus 

Let C be any d-variate copula and let Ui,U2, ■ ■ ■ be an i.i.d. sequence of random vectors 
with common cumulative distribution function C. Let q;„ be the multivariate empirical 
process in equation (2.1). Consider the oscillation modulus defined by 

Af„(a) = sup{|q;„(u) — a„(u)|: w,ue[0,l]'', |uj — | < a for all j} (A.l) 

for a € [0, oo). Define the function tjj : [— 1, oo) — > (0, oo) by 

il;{x)= 2x-^{{l + x)log{l+x)~x}, a;e (-l,0)U(0,oo), (A.2) 

together with 1) = 2 and ip{0) = 1. Note that ip is decreasing and continuous. 

Proposition A.l (John H. J. Einmahl, Hideatsu Tsukahara). Let C be any d- 

variate copula. There exist constants Ki and K2, depending only on d, such that 

Pr{Af„(a) > A} < ^ exp(^:^^ ( ^) ] 
for all ae{{),l /2] and all A £ [0, 00) . 

Proof. In Einmahl [7], Inequality 5.3, page 73, the same bound is proved in case C is 
the independence copula and for a > such that 1/a is integer. As noted by Tsukahara, 
in a private communication, the only property of the joint distribution that is used in the 
proof is that the margins be uniform on the interval (0, 1): Inequality 2.5 in Einmahl [7], 
page 12, holds for any distribution on the unit hypcrcubc and equation (5.19) on page 72 
only involves the margins. As a consequence, Inequality 5.3 in Einmahl [7] continues to 
hold for any copula C. Moreover, the assumption that 1/a be integer is easy to get rid 
of. □ 
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